Observation of dynamical instability for a Bose-Einstein 
condensate in a moving ID optical lattice 
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We have experimentally studied the unstable dynamics of a harmonically trapped Bose-Einstein 
condensate loaded into a ID moving optical lattice. The lifetime of the condensate in such a potential 
exhibits a dramatic dependence on the quasimomentum state. This is unambiguously attributed to 
the onset of dynamical instability, after a comparison with the predictions of the Gross-Pitaevskii 
theory. Deeply in the unstable region we observe the rapid appearance of complex structures in the 
atomic density profile, as a consequence of the condensate phase uniformity breakdown. 
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Ultracold atoms in optical lattices have proven to be a 
rich field of investigation for both fundamental and ap- 
plicative issues of quantum mechanics. Periodic poten- 
tials are very well known in solid state physics, where 
the basic description of the system is the Bloch theory 
for a gas of noninter acting particles, such as electrons 
in crystalline solids Cold atoms in optical lattices 
have opened the possibility to investigate effects not pre- 
viously observable on ordinary matter crystals, such as 
Bloch oscillations, Wannier- Stark ladders and Landau- 
Zener tunnelling |2|. Bose-Einstein condensates (BECs) 
are particularly well suited for the investigation of these 
phenomena, due to the small momentum spread and the 
large coherence length. However, in a high density sam- 
ple, such as a trapped BEC, nonlinearities induced by 
interactions among the atoms can complicate the pic- 
ture and lead to a number of new effects, such as the 
observation of new quantum phases |3(, the generation 
of solitonic propagation 3 and the observation of differ- 
ent kinds of instabilities. In the recent years many pa- 
pers have investigated the latter topic, both theoretically 
U H H S H Gil and experimentally [H E O . Many 
of these works have suggested that the BEC superflow 
can be broken not only by the mechanism of energetic 
(Landau) instability, but also by dynamical instability, 
occurring for condensates with repulsive interactions in 
periodic potentials. 

In this Letter we report on the unambiguous obser- 
vation of dynamical instability for a trapped BEC in a 
moving ID optical lattice. We have measured the char- 
acteristic rates for this instability to occur, mapping the 
quasimomentum space for both the lowest and excited 
energy bands. This quantitative analysis gives important 
information on the timescales limiting the investigation 
of coherence effects in a BEC moving in a periodic po- 
tential. 

Dynamical instability is a peculiar feature of nonlinear 
systems. It occurs when the eigenspectrum of the excita- 
tions of the system exhibits complex frequencies. In this 
case, arbitrary small perturbations of the wavefunction 
may grow exponentially, eventually leading to the de- 
struction of the initial state. The conditions for this kind 



of instability are satisfied by a BEC with repulsive inter- 
actions in a periodic potential. When the height of the 
periodic potential is not too large, the superfluid ground 
state of a BEC in an optical lattice is well described by 
the Gross-Pitaevskii equation (GPE) 
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where \I/ is the complex BEC order parameter, Vtrap(r) 
is the harmonic trapping potential, s is the height of the 
periodic potential in units of the recoil energy Eji = 
h 2 k 2 /2m, n/k is the periodicity of the lattice in real 
space, and a is the scattering length. Among the station- 
ary solutions of Eq. there are the usual Bloch waves, 
i.e. plane waves with quasimomentum q and band index 
n, modulated in space by functions having the same pe- 
riodicity of the lattice. A linear stability analysis of the 
GPE shows that, in certain regions of the quasimomen- 
tum space, these solutions are not dynamically stable 
and an exponential growth of perturbations may start 
IE IE EH • The onset of dynamical instability has been 
considered as a possible explanation for the disruption of 
the atomic superflow in a past experiment [l2|, concern- 
ing large amplitude dipole oscillations in the magnetic 
trap + static optical lattice. In this work we quanti- 
tatively investigate the unstable dynamics of a trapped 
BEC in a moving optical lattice in the regime of low lat- 
tice heights. The novelty of this work, with respect to the 
previous experiments [Tj, [13 , is that the implementation 
of a lattice moving at constant velocity allows us to ac- 
curately set the quasimomentum of the condensate and 
to make a precise investigation of the stability regimes in 
the full Brillouin zone and for different energy bands. 

The experiment is performed with an elongated 87 Rb 
Bose-Einstein condensate produced in a Ioffe-Pritchard 
magnetic trap by means of forced RF evaporation. The 
trapping frequencies are uo z = 2tt x 8.8 Hz axially and 
lu± = 2tt x 91 Hz radially, with the axis of the trap ori- 
ented horizontally. Our typical condensates are made 
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FIG. 1: Loss rates for a trapped BEC loaded in a moving 
optical lattice with s = 0.2. The vertical line corresponds 
to the calculated threshold for the onset of dynamical insta- 
bility 10]. The images show the density distribution of the 
expanded cloud. Near the zone boundary, where instability is 
faster, we observe the appearance of some complex structures, 
evidencing the loss of coherence in the BEC. 



of ~ 3 x 10 5 atoms in the hyperfine ground state \F = 
1; mp = —1 >, with a peak density n ~ 1.2 x 10 14 cm -3 . 
The optical lattice is created by the interference of two 
counterpropagating laser beams derived from the same 
Ti:Sa laser operating at A = 820 nm, far detuned with 
respect to the Rb Dl line at A = 795 nm. The lattice 
beams are aligned along the symmetry axis of the con- 
densate and are only slightly focused (400 fim diameter), 
so that the optical radial confinement may be completely 
neglected. The frequencies of the two beams are con- 
trolled by two acousto-optic modulators (AOMs) driven 
by two phase-locked radiofrequency generators in order 
to provide a stable detuning 5v between the two beams. 
The resulting interference pattern, averaged over the op- 
tical frequencies, is a standing wave moving at velocity 
v = (\/2)5v. Once the condensate has been produced, 
we adiabatically switch on the moving lattice loading the 
condensate in a state of well defined quasimomentum 
q = mv/h and band index n We let the BEC evolve 
in this potential for a variable time, then we switch off 
both the magnetic trap and the optical lattice and, after 
an expansion of 28 ms, we ima ge t he atomic cloud along 
the radial horizontal direction |l5j . 

The number of atoms remaining in the condensate de- 
creases exponentially as a function of the time spent in 
the periodic potential. We have measured the lifetime of 
the condensate for different values of the quasimomen- 
tum q and for different energy bands. Since even a tiny 
thermal component may seriously limit the BEC lifetime, 
we have used an RF-shield in order to remove the hottest 
atoms produced by heating of the atomic sample. In this 
way we measure lifetimes of the order of ~ 10 s in the lat- 



tice at v — 0, with no discernable thermal fraction even 
on long timescales. 

The experimental results are analyzed according to the 
GPE shown in Eq. (1). As a matter of fact, one can 
avoid the complication of the full 3D theory by using the 
Non Polynomial Schrodinger equation (NPSE), a simpli- 
fied ID model that includes an effective radial-to-axial 
coupling The NPSE provides a more realistic de- 

scription of the actual system with respect to the simple 
one dimensional GPE used previously p , yielding an es- 
timate of the instability thresholds and of the growth 
rates of the most unstable modes in nice agreement with 
that of Eq. (1) in the regime of lattice heights considered 
here (see 10] for details). 

In Fig. n we show the measured loss rates (inverse of 
the lifetime) for the BEC in an optical lattice with s = 0.2 
as a function of q. With increasing q, from the bottom of 
the first band to the zone boundary, the lifetime changes 
dramatically, spanning three orders of magnitude from 
« 10 s to « 10 ms. In particular, we observe a disconti- 
nuity around q = 0.55g# (where qs = 27r/A is the bound- 
ary of the first Brillouin zone), in good correspondence 
with the calculated threshold for the onset of dynami- 
cal instability (vertical line). Deeply in the dynamically 
unstable regime, as shown in the images of Fig. ^ we 
observe the appearance of some complex structures in 
the expanded BEC density profile, suggesting fragmenta- 
tion of the Bloch wave. These interference-like structures 
are particularly visible for higher lattice heights and near 
the zone boundaries, where the evolution of instability is 
faster. In Fig. we present two sequences of images 
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FIG. 2: A) Evolution of the expanded atomic density profile 
for a lattice height s = 1.15 and two different quasimomentum 
values q = OAqs and q = 0.55^ below and above the thresh- 
old for dynamical instability. Note the different timescales. 
B) Evolution of the expanded atomic density profile for a 
variable time spent in the pure magnetic trap, following 5 ms 
spent in an optical lattice with s — 1.15 and q = 1.30#b. In 
all these pictures the lattice is directed from top to bottom. 
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FIG. 3: A) Experimental loss rates for a BEC loaded into 
the first three energy bands of a moving optical lattice with 
s = 0.2. B) Theoretical growth rates of the most dynamically 
unstable modes obtained from a linear stability analysis of 
theNPSE[T3. 



showing the time evolution of the atomic density profile 
for a lattice height s = 1.15 and two different values of 
quasimomentum below and above the threshold for dy- 
namical instability. This behavior cannot be attributed 
to a simple heating of the sample. On the contrary, it 
may reflect the creation of phase domains induced by the 
growth of instabilities, which break the phase uniformity 
of the BEC p. When these structures are not just weak 
perturbations of the density profile, it is difficult to mea- 
sure the number of atoms in the condensate directly. In 
order to measure it more precisely, we let the cloud evolve 
in the pure magnetic trap, allowing relaxation of the ex- 
citations. We have measured typical relaxation times of 
the order of « 500 ms, after which the BEC recovers its 
smooth density profile, as shown in Fig. [2f>. This obser- 
vation suggests the existence of mechanisms that, once 
the cause of instability is removed, allow the system to 
come back to the ground state, damping the excitations 
and restoring the initial coherence. 

Our experimental technique allows us to study the dy- 
namics in the higher energy bands and compare the pe- 
culiar behavior in each band with the results of the the- 
oretical model. The remarkable agreement between the 
theory and the experiment enforces that the observed 
phenomena are indeed due to dynamical instability. In 
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FIG. 4: A) Experimental loss rates for a BEC loaded into the 
first three energy bands of a moving optical lattice with s = 
1.15. B) Theoretical growth rates of the most dynamically 
unstable modes obtained from a linear stability analysis of 
theNPSE[l(J. 



Fig. OlA we show the experimental loss rates measured 
in the lowest three energy bands as a function of the 
quasimomentum q for s = 0.2. These rates are com- 
pared in Fig. [3f> with the theoretical growth rates of the 
most dynamically unstable modes, obtained from the lin- 
ear stability analysis of the NPSE |l£|. Increasing the 
height of the optical lattice from s = 0.2 to s = 1.15, as 
shown in Fig. the picture in the higher bands starts 
to develop asymmetric features around the zone edges. 
This is particularly evident crossing the boundary be- 
tween the second and the third Brillouin zone, where we 
observe a marked change in the instability rates. This 
nontrivial behavior cannot be attributed to single parti- 
cle band structure considerations (for example, residual 
nonadiabaticites or interband transitions). Indeed, com- 
paring this feature with the theoretical calculations for 
the growth rates of the most unstable modes, we do ob- 
serve the same distinctive shape. The experimental loss 
rates should not be quantitatively compared with the the- 
oretical growth rates, since these two quantities have dif- 
ferent physical meanings: the first measures how fast the 
atoms are removed from the condensate; the second is the 
rate at which the most unstable mode grows in the linear 
regime, i.e. at the onset of dynamical instability. Out of 
the linear regime, when excitations have grown and the 
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unstable modes are not just weak perturbations of the 
carrier Bloch wave, the dynamical evolution of the sys- 
tem cannot be explained with this perturbative approach 
and the full solution of the time-dependent 3D-GPE is re- 
quired. However, the remarkable similarity between the 
experimental and theoretical curves indicates that the 
onset of the instability produces a significant imprinting 
on the subsequent dynamics of the system. 

In the experiment, the BEC is trapped in the harmonic 
magnetic potential, in order to have a high density sam- 
ple and long observation times. On the other hand, the 
presence of the harmonic potential may induce a varia- 
tion of the quasimomentum in time, thus changing the 
conditions for instability. In a semiclassical approach, a 
force F acting on the system results in a variation of the 
quasimomentum q according to the law F = dq/dt. In- 
tegrating this equation of motion with the velocity spec- 
trum v(q) of the Bloch bands, one can demonstrate that 
the quasimomentum q makes small oscillations at the 
axial trap frequency (rescaled with the effective mass) 
[~L7^ . We have taken this effect into account by including 
horizontal error bars in Figs. |3| and 21 representing the 
maximum variation of q during the experiment. In real 
space these oscillations correspond to a micro-motion of 
the BEC around the center of the trap with an ampli- 
tude much less than 1 /am. Actually, this effect becomes 
important only when the Bloch bands significantly differ 
from the free particle energy, i.e. near the zone bound- 
aries and for higher lattice heights. 

Repeating the experiment in the presence of a small 
thermal component we observed that, even for small val- 
ues of quasimomentum (q > 0.05#b), the atomic sample 
is completely destroyed in a time much shorter than the 
lifetime measured for the pure BEC. We attribute this 



behavior to the onset of energetic instability, that occurs 
in the presence of dissipative processes 5, 101. The resid- 
ual thermal fraction surrounding the BEC can provide a 
mechanism for such a dissipation, thus triggering the ac- 
tivation of energetic instability. The study of this regime 
will be the subject of further experimental work. 

In conclusion, we report the observation of dynami- 
cal instability for a harmonically trapped Bose-Einstein 
condensate in the presence of a ID optical lattice. The 
ability to precisely control and change the lattice veloc- 
ity has allowed us to carry out a detailed investigation of 
the instability regimes in the full Brillouin zone for dif- 
ferent energy bands. We have quantitatively studied the 
lifetime of the condensate in such a potential, finding a 
good agreement between the experimental loss rates and 
the theoretical rates for the growth of excitations at the 
onset of dynamical instability. Deeply in the dynami- 
cally unstable regime we have observed the appearance 
of complex structures in the expanded density profile, a 
signature of a loss of coherence in the atomic sample. 
These observations clearly identify dynamical instability 
as the main mechanism responsible for the disruption of 
the matter wave superflow in the regime of small lat- 
tice heights. The quantitative analysis presented in this 
letter gives important information on the characteristic 
timescales which limit the investigation of the coherence 
properties of a BEC moving in a periodic potential. 
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